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Abstract. The role of shake-off for double ionization of atoms by a single photon 
with finite energy has become the subject of debate. In this letter, we attempt to clarify 
the meaning of shake-off at low photon energies by comparing different formulations 
appearing in the literature and by suggesting a working definition. Moreover, we 
elaborate on the foundation and justification of a mixed quantum-classical ansatz for 
the calculation of single-photon double ionization. 
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It is well known that the double-to-single cross section ratio for ionization of atoms 
by a single photon does not vanish at high photon energies. Rather, it approaches a 
finite constant which can be explained in the framework of a sudden approximation by 
the so-called shake-off mechanism. While shake-off is well defined in the asymptotic 
high-energy limit, its meaning at finite energies is less clear and has been the subject 
of debate recently. In fact, a number of definitions can be found in the literature, 
e.g., [H 12 El lU El El • Some of these definitions are based on formal diagrammatic 
perturbation expansion techniques [SI E], others on more general physical arguments 
[211], others on a simple extension of the sudden approximation idea to finite energies 
PU El d- Most of them (with the exception of [3j, see below) have in common that 
they approach the well-known asymptotic expression at high energies. At low and 
intermediate energies, however, they may differ markedly (e.g., some show a monotonic 
dependence on energy while others have a maximum at some finite energy, some even 
exceed the total double-to-single cross section ratio measured experimentally). Thus, so 
far no unique definition for shake-off at finite energies exists, and it is not obvious what 
"the best" definition might be. On the other hand, in particular in connection with 
the interpretation of experimental data, attention has been given to the question what 
physical mechanisms dominate double ionization and at which energies the different 
mechanisms are important [S]. 

One may argue that a satisfying definition of shake-off would be one based on 
physical principles in addition to mathematical rigor. Hence, a "good" definition should 
separate shake-off as much as possible from other ionization mechanisms. Clearly, such 
a separation is not strictly possible in the presence of other available routes to ionization 
and can only be approximate, which makes the discussion of a shake-off mechanism a 
somewhat delicate issue. In comparing calculations with experiments, e.g., one should 
always keep in mind that no strict one-to-one correspondence between a shake-off 
mechanism as an approximate physical picture and a separate calculation of shake- 
off can be expected due to the neglect of interference between possible decay routes. 
Nevertheless, since such simple physical pictures, to the extent of their applicability, 
can be very valuable for our intuitive understanding of physical processes, a definition 
separating shake-off from "non- shake-off " would seem most rewarding conceptually. One 
such definition has recently been given by Schneider et al. [Z] (hereafter referred to as 
SCR), where the single-photon double ionization process has been described in terms of 
two separate contributions, namely "shake-off" and "knock-out". The method used in 
SCR was shown to lead to excellent agreement with experiment and ab-initio calculations 
for double ionization from the ground state jZUHl, and very recently also from excited 
states [rU], of helium. Thus, it suggests itself as a more or less natural "operational" 
definition of shake-off in the framework of the "half-collision" picture of single-photon 
multiple ionization HJ . 

The calculation reported in SCR starts from a mixed quantum-classical ansatz 
that is based on the separation of the photoabsorption process (which is not treated 
explicitly in the calculation) from the subsequent evolution of the system. It treats this 
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evolution (i.e. the redistribution of the energy between the two electrons) in the spirit of 
an (e,2e)-like process with the additional possibility of shake-off. Such a "half-collision" 
picture has been originally suggested by Samson [TT] and elaborated by Pattard and 
Burgdorfer j^, allowing for shake-off processes which are not taken into account in 
Samson's original model. In the SCR ansatz, the (e,2e)-like ("knock-out") part of the 
cross section is calculated using a classical trajectory Monte Carlo method, to which 
the shake-off as a purely quantum mechanical process is added on top. In this spirit, 
shake-off is introduced as a more or less ad-hoc quantum correction to an essentially 
classical treatment. Here, we start from a fully quantum mechanical expression and see 
which kind of approximations lead to an SCR-like ansatz. In this way, further insight 
into the validity of the ansatz, concerning both technical details of the calculation as 
well as the approximate separation of physical mechanisms (shake-off and knock-out), 
can be obtained. 

In ref. a Born series for the transition amplitude from the ground state ipi 
to a final state xpf of a two-electron target following single-photon absorption has 
been derived. It was shown that, under the assumption of negligible electron-electron 
correlation in the ionized final state, the transition amplitude can be written as 

a/i = -2m S {E f -E t -u) (^f| ^1 - i Jolt e iHot T ee e" iHot j V pe \^) . (1) 

In the above equation, V pe is the photon-electron interaction, usually taken in dipole 
approximation, Hq is the final-state Hamiltonian Hq = H at — V ee and T ee denotes the 
Coulomb T-matrix for electron-electron scattering, ip^p is an eigenfunction of Hq, i.e. a 
product of two one-electron states, due to the assumption of vanishing electron-electron 
correlation in the final state (where at least one electron is ionized), while ipi is the 
fully correlated initial (ground) state of the target. Introduction of a complete set of 
intermediate states then allows for a separation (on the amplitude level!) of the initial 
photon absorption from the subsequent propagation 

a fl = -2m 5(E f -Ei-u) ^^/V+IVO (^a\V pe \^) , (2) 

a 

where the notation 

S+ = l-i f dt e lHot T ee e~ iHot (3) 



is motivated by its resemblance of a conventional scattering S-matrix. Note, however, 
that S+ is not strictly an ^-matrix for electron-electron scattering since the time 
integral in equation Q is restricted to positive t, i.e. S + corresponds to a /ia//-collision. 
Furthermore, let us choose the complete set {ip a } in such a way that 

W1,2)e t JM= (4) 
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is contained in this set. From the orthogonality condition for the basis states, all other 
basis states are then orthogonal to V pe ipi. Thus the sum over intermediate states in 
equation (J2J collapses to a single term and we can write 

a fi = -2m 5(E f - Ei-u) (V>/ \S + \il) abs ) {ip a b s \V pe \ipi) . (5) 

The photon absorption probability is then given by a sum over all final states ipj of 
the transition probability per unit time into the state ip^p 

P ahs = 2 7 rJ2 5 (Ef-E i - oo) IS+IVV*,)! 2 \(^abs\V pe \^)\ 2 . (6) 

/ 

On the other hand, it is also directly given by 

P abs = 2irJ2HE f -E l ~uj) \(tfj f \V pe \^)\ 2 . (7) 
/ 

(Note that the ipf in equation (J7j) are eigenfunctions of the full atomic Hamiltonian H at 
including electron-electron interaction, in contrast to the final states in equation (JHJ).) 
From equation (J3J), however, it immediately follows that 

|(^& s |^#i>| 2 = (V pe ^i\V p ^i) = J2\&f\ v vM\ 2 > ( 8 ) 

/ 

which in general does not coincide with the expression (JZj) involving an additional delta 
function. Hence, it can be seen that it is precisely the off-shell (i.e. off the final- 
state energy shell) part of ip a b s which prohibits an exact factorization of the transition 
probability into a photon absorption probability and an "energy redistribution" part. 
Note that this can also be seen from equation (j2j) directly if the set of intermediate 
states tp a is chosen as eigenstates of H . Then 

oo 

(^f = S fa -tJdt e « E f- E ** ^f\T ee \^ a ) (9) 



and from 

oo 

/ dt e^-^X = nd (E f - E a ) + i (10) 



it becomes clear that the off-shell part of ip a b s is a consequence of time ordering 
i.e. that requiring the photon to be absorbed first restricts the time integral in (J3J) to 
positive t. 

For the remainder of this discussion, let us neglect the photon absorption process 
and focus on the second step of the ionization process, namely the "half-collision" part 
of equation (j5J) 



a Labs = V2^5(E f -Ei-u) (tpf\l -ijdt e lHot T ee e"^!^) • (11) 

o 

(The splitting of the factor 2ti is motivated by the fact that the resulting shake probability 
to be discussed below reduces to the correct asymptotic form at high energies.) af >abs 
is seen to consist of two parts, namely the interaction free unity operator "1" and 
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the operator "T" involving electron-electron interaction. Naturally, the former can be 
associated with a shake process while the latter corresponds to a "knock-on" (we use the 
expressions shake and knock-on for any final state and the terms shake-off and knock-out 
for doubly ionized states as in SCR). Hence, we propose 



a 



S = ./^ KIT?. _ T? ^ /„/,(°) 

f,abs _ 



27rS(E f -E t -uj)(ijy ) \^ abs } (12) 



as a working definition for the shake amplitude at a finite excess energy E = E{ + u. 
However, the shake and knock-on contributions are summed on the amplitude level. To 
arrive at the SCR ansatz, the additional approximation of an incoherent summation of 
shake and knock-on has to be made:):. The error introduced by this approximation is 
at most of the order of the smaller of the two contributions, i.e. it goes to zero in the 
high- as well as low-energy limit and could only contribute significantly at intermediate 
energies. Even there it was found in SCR that the error is of the order of a few percent 
only (at least for the double-to-single ionization ratio of helium). One would speculate 
[3J that this is to a large extent due to the population of different final states by the 
two mechanisms. For shake, e.g., the "shaken" electron will be in an s-state, while the 
knock-on mechanism will also populate higher angular momentum states. Calculations 
of angular-differential cross sections should shed further light on this question. 

From equation ()12j) . the probability for a shake process to a final state ipf per unit 
time is found to be 

Pl bs = KE S -E) |<^°W| 2 . (13) 

With the definition of tp a bs, equation (J3J), this is more explicitly written in terms of the 
initial state as 



|(V> ( f 0) |K e V>i)| 2 

p '^ i(E '- E) f^k- (14) 



This expression differs somewhat from the one given by Aberg pQ (equation (7) of 
SCR). In contrast to the former, it contains the "photoabsorption operator" V pe . It 
seems that the current definition equation ()14j) is preferable since it arises naturally 
from the preceding arguments: The sudden approximation underlying the shake-off 
picture is with respect to the electron-electron interaction, not with respect to the 
photon absorption. As noted by SCR, if the photoelectron is in an s-state initially, it 
will be in a p-state after absorption of the photon. That is, to the extent that the 
dipole approximation is valid for the photon-electron interaction and that the single- 
particle angular momentum I is a good quantum number, the ip v defined by Aberg 
is identically zero. It should be noted, however, that both expressions lead to the 
same high-energy limit. In this limit, as well as (ipf^ \V pe ipi) both become 

proportional to (4> ef (r 1 )\ip i (r 1 , r 2 = 0)) (where cf> denotes the one-electron state of the 
shaken electron and e its energy) with different prefactors. Since they appear equally 

| In addition to that, the knock-on part has been obtained from a classical CTMC calculation in SCR. 
While such a treatment is frequently employed in the study of atomic collision processes, an evaluation 
of its quality is beyond the scope of the present Letter. 
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in the numerator and the denominator of (|14|) they cancel out, leading to the same 
high-energy limit. The same is true for all energies if, as in SCR, product wavefunctions 
are used for the initial state, or if the PEAK approximation is employed [UIHE], i-e. if 
the absorption of the photon is assumed to happen always at the nucleus at any excess 
energy E. In this case one arrives at the "natural" definition (equation (8) of SCR) 

P* abs = 9 (E - e f ) m f m 2 (15) 

(where 9 is the unit step function), i.e. the overlap of two one-electron wavefunctions. (In 
the case of using the PEAK approximation for a correlated initial state <f>i = V^( r i> r 2 = 
0)/(^(r 1 ,r 2 = 0)|^(r 1 ,r 2 = 0)) 1 /2.) 

As argued above, the successful application of the SCR method showing excellent 
agreement with experiment and ab-initio calculations suggests the adoption of equation 
f!14|) as a good "operational" definition of shake-off. In this spirit, shake-off may be 
phrased vaguely as the part of the double ionization that is absent in a full collision (due 
to the orthogonality of initial and final states), or, more precisely, that part which does 
not involve an electron-electron interaction explicitly (of course, interaction is implicit in 
the correlation present in the initial state, without which there would be no shake-off). 
It should be noted that the quality of this definition depends on the observation that 
there is very little interference between shake-off and knock-out, as discussed above. 
The fact that this is not strictly true leads to some problems e.g. at very low energies, 
where this separation would lead to a linear dependence of the double ionization cross 
section on the excess energy, in contrast to the well-known Wannier threshold power 
law [IB]. 

On the other hand, one might want to adopt a maybe more "physical" definition 
of shake-off on the basis of the intuitive picture of a time-dependence of the effective 
one-electron potential the "shaken" electron feels. For asymptotically high energies, it 
is the sudden change of this potential that leads to a relaxation of the electron which 
is not in an eigenstate anymore after the potential has changed. In this sense, the 
change in the effective potential does not occur suddenly anymore at finite energies, but 
rather over a timescale given by the velocity of the outgoing photoelectron. This is the 
basis for the definition of shake-off adopted in [2j , where an expression has been derived 
from time-dependent perturbation theory, and also the rationale behind the somewhat 
ad-hoc formulation used in [I] which is motivated by a Rosen-Zener-like expression 
for diabatic transition probabilities familiar from ion-atom collisions. It should be 
noted that, in contrast to ([Tlj) . both of these expressions show an exponential decrease 
towards threshold, so that the Wannier threshold law is recovered. While in the SCR- 
expression the probability to be shaken into a specific final state does not depend on 
the rate of change of the potential (i.e. the velocity v of the photoelectron) as long as 
it is energetically allowed, this is different in [21 E] where these probabilities depend 
exponentially on v 2 and v , respectively. In view of these differences, it is surprising to 
see that the numerical values resulting from these different definitions are in fact rather 
similar, as is demonstrated in figure 1. This leads us to suggest f!14|) as a good working 
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Figure 1. Comparison of different shake-off expressions for the case of the helium 
ground state. Solid line: Schneider et al. 7 , dashed: Thomas dot-dashed: Pattard 
and Burgdorfer 0|; the thick solid line shows the asymptotic E — > oo limit. 

definition of shake-off at finite energies. It agrees qualitatively with other natural and 
maybe more "physical" definitions and, moreover, it has been shown to lead to a very 
good approximate separation of mechanisms into shake-off and knock-out. It is easy to 
calculate and, as compared to [21 H] , has the advantage that it does not contain any free 
parameters (such as a characteristic range of interaction as in [2] or an effective impact 
parameter as in However, one should always keep in mind that the significance of 
(fTIJ) as an independent physically meaningful quantity is limited, as discussed, e.g., in 
connection with the behaviour of the cross section near threshold. 

At this point, a short comment on a comparison with other available definitions 
of shake-off seems to be in place. So far we have given a comparison with [21 E] only, 
and used the qualitative agreement with the values calculated from the expressions given 
there as an argument in favour of the current definition. However, as argued before, some 
other definitions found in the literature lead to shake-off values rather different from 
the present. The shake-off calculated from many-body perturbation theory (MBPT) 
[3j, e.g., is found to have a completely different behaviour. Its shape (as a function 
of energy) as well as the high-energy limit reached, even the answer to the question 
whether it is the dominant process at high energies or not, depend on the choice of the 
gauge in which the corresponding diagram is calculated. This is not too surprising in 
view of the fact that only the sum of all first-order diagrams (shake-off plus ground-state 
correlation plus two-step one) has a well-defined and gauge invariant meaning. Hence, 
the meaning of shake-off is well defined within MBPT, however it is not claimed to have 
any independent physical meaning of its own. In this sense, it is not a helpful quantity if 
one wishes to discuss approximate physical mechanisms. Another definition, originally 
formulated by Aberg pQ, has recently been used by Shi and Lin [H] to calculate shake-off 
double ionization of the helium ground state. Their result for the double-to-single cross 
section ratio is found to be significantly larger than the latest experimental data for 
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the total ratio, i.e. including all possible decay routes (e.g. shake-off and knock-out). 
From jH] it is not entirely clear how much of this "overshooting" has to be attributed 
to a poor choice of the ground state wavefunction (leading to an asymptotic high- 
energy limit which is somewhat too large) and how much would still be observed using 
a more accurate initial state. Assuming that this effect persists it would be obvious 
that again no physical meaning can be ascribed to this definition of shake-off, since 
shake-off alone would already be larger than the sum of all mechanisms. In any case, a 
further discussion would have to await a corresponding formulation of "non- shake-off " , 
since it is only the sum of all possible ionization mechanisms which can directly be 
compared with experiment. Finally, Kheifets has proposed a definition of shake-off 
where the diagonal part of the T-matrix contribution to the convergent close-coupling 
model is absorbed into shake-off. In his calculations for helium, it was found that with 
this definition of shake-off the total cross section ratio approaches the shake-off value 
quickly, and non-shake-off becomes negligible at about 100 eV excess energy. However, 
the last panel of figure 2 of [Sj shows that for lower energies shake-off alone again exceeds 
the total ratio. Hence, once more one has to conclude that the meaning of shake-off as 
defined in as an independent physical mechanism is limited. 

In summary, we have argued that no unique definition for shake-off at finite energies 
exists. Nevertheless, we propose equation (fPIj) as a good "operational" definition. 
Clearly, when shake-off at finite energies is discussed (in particular in the sense of 
an approximate physical mechanism in connection with experiments), care has to be 
taken of the precise meaning of the term, i.e. its actual definition adopted in each case. 
In addition to our discussion of shake-off, we have indicated a way towards a rigorous 
derivation of the SCR ansatz for calculation of double ionization by relating it to a 
perturbation expansion starting from a full quantum mechanical point of view. 

T.P. would like to thank Joachim Burgdorfer for stimulating his interest in the 
problem of single-photon multiple ionization, as well as many very helpful discussions. 
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